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Coulomb interaction and instability of CE-type 
structure in half-doped manganites 

In their Letter, [1] den Brink, KhaliuUin, and Khom- 
skii proposed theoretically that the one-dimensional fer- 
romagnetic zigzag chains in CE phase in half-doped man- 
ganites play an essential role in forming the orbital order- 
ing, and, more surprisingly, the on-site Coulomb interac- 
tion U between electrons with different orbitals leads to 
experimentally observed charge ordering. In this Com- 
ment, I point out that the strong U will destroy the sta- 
bility of CE-type phase, which is stable in a very narrow 
regime in the parameter space for electronic model. 

The same notations as in Ref. [1] are employed here. 
In Ref. [1] , the case otU = was analyzed and CE phase 
was found to be stable when the superexchange coupling 
J > 0.1524t (t is the hopping integral). In the CE phase 
the one-dimensional zigzag chain consists of two types 
of lattice sites, the bridge and corner sites. The bridge 
site is occupied singly by electrons with certain orbitals 
while the corner sites can be occupied doubly. With- 
out U, the charge distribution along the zigzag chain is 
uniform, i.e., 1/2. Finite U will push the charge away 
from the corner sites to the bridge sites to form charge 
ordering. The theory explains the patterns of three or- 
derings almost perfectly. We observe that the U will also 
increase the energy of the CE phase. The rod-type (i.e., 
C-type) phase consists of one-dimensional ferromagnetic 
straight chains. The chains are occupied at most singly 
on each site. Thus the U does not change the energy 
of the C-type phase. This fact leads to the instability 
of CE phase with respect to C-type phase for finite U. 
We evaluate the electronic energy for CE phase by us- 
ing exact diagonalization method. The energies from the 
magnetic coupling are the same for both C-type and CE 
phase. At [/ = 0, E^^ = -0.6953t for CE-type phase 
and i?^„ = — 0.6366i for C-type phase. However, at 
U = +0O, E^-^ = — 0.61t and the energy of C-type phase 
keeps unchanged. In other words, the CE phase is unsta- 
ble with respect to C-type phase when electron correla- 
tion becomes stronger. Numerical results for energies of 
CE phase via the U is shown in Fig. 1(a). The calculation 
is performed on 8-, 12- and 16-site chains, respectively. 
The size-effect is very small. The main reason is that 
the half-doped zigzag chain is a band insulator. This is 
also seen from the calculation of E^J^ for a finite chain 
when U = 0. Its value for 16 sites is -0.695328t, which 
is equal to the value of an infinite chain up to the order 
of 10~^.[2] At about U ~ 5t, the CE-type phase has a 
higher energy than C-type phase. To establish a compre- 
hensive phase diagram, we have to consider other possible 
phases. Except for conventional anti- and ferromagnetic 
phases, a novel type of phase is taken into account. The 
phase consists of two-site ferromagnetic valence bonds 
(VB), and the bonds couple antiferromagnetically. Elec- 
trons are localized within the bonds. The electronic part 
of the state in a valence bond is :^(cl,a,o- + ^],a.a) 
and the local spins are parallel to spin of electron. The 



energy per site is E^^^ + E^^^ = —0.5t — J, which is in- 
dependent of U (we take the magnetic energy as zero 
for CE phase). The ferromagnetic phase has energy 
—t + 2J for three dimension and — 0.919i -I- 2J. Thus 
when [/ = and 0.152t < J < 0.195t for three dimension 
and 0.117i < J < 0.195t for two dimension, CE phase 
has a lower energy. The phase diagrams for J7 = and 
U — +00 as shown in Fig. 1(b) and (c) are evaluated on a 
4x4 lattice with a periodic boundary condition by means 
of combination of mean field theory and exact diagonal- 
ization. We assume the magnetic structure and claculate 
the energy by exact diagonalization. FM and C-AF are 
assumed to contain two sublattices and evolve into G-AF 
when spins on two sublattices are not parallel for large J. 
Thus, the CE phase does not appear in the ground state 
for the model in Ref. [1] in the case of strong correlation. 

In doped manganites, U is usually estimated to be 
much larger than t. In Ref. [1], U k, lOt. The on-site U 
alone can produce the experimentally observed pattern 
of charge ordering, but cannot stabilize the CE phase. 
To solve this issue finally, we have to take into account 
other interactions, such as the long-range Coulomb inter- 
action, Jahn- Teller distortion, and physics of topological 
berry phase. For example, the effect of finite large 
leads to an attractive particle-hole interaction, which fa- 
vors to stabilize the charge ordering. [3] 
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FIG. 1. (a) The electronic energy of CE phase per site via 
the on-site interaction; (b) and (c) the phase diagrams on a 
square lattice for (7 = and -foo, respectively . 

Shun-Qing Shen 

Department of Physics, The University of Hong Kong, 
Hong Kong and Institute of Physics, Chinese Academy 
of Science, P. O. Box 603, Beijing 100080, China 



1 



^ J. van den Brink, G. KhaliuUin, and D. Khomskii, Phys. 

Rev. Lett. 83, 5118 (1999). 
^ The finite-size effect for E^^^ is relatively large. i5£„(12) = 

-0.6220f while £;g„(+oo) = -0.6366f. So the calculated 

value Uc from Fig. 1(a) should be more reasonable than 

that from 12 sites. 
^ S. Q. Shen and Z. D. Wang, Phys. Rev. B 59, 14484 (1999); 

61, 9532 (2000). 



2 



